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Abstract: We compute the complete set of Feynman Rules producing the Rational Terms 
of kind R2 needed to perform any QCD 1-loop calculation. We also explicitly check that in 
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1. Introduction 

In the last few years, the problem of computing one-loop amplitudes efficiently has been 
attacked by several groups. Standard techniques, such as the Passarino-Veltman tensor 
reduction and its many variances Q-Q, have been used for many years, and have produced 
a great deal of useful results [||. Nowadays new developments, in which the amplitude is 
directly reconstructed are widely used. Such approaches rely on the fact that the basis of 
one- loop scalar integrals is known in terms of Boxes, Triangles, Bubbles and (in massive 
theories) Tadpoles, so that any one-loop amplitude A4 can then be written as: 

M = d i Box * + Ci Trian § le i + bi Bu bblei + ^2 a i Tadpolej + R , (1.1) 

i i i i 

where dj, Cj, bi and en are the coefficients to be determined and R is a remaining piece, 
called Rational Term (RT). The first attempt in this direction was the Unitarity approach 
Q, by Bern, Dixon and Kososwer, in which two-particle cuts are performed on the one-loop 
amplitude (or, equivalently, two tree-level amplitudes are glued) in order to get information 



on the coefficients in eq. 1.1. The method produced many useful - mainly analytical 



results, especially for massless theories M, but a systematic way to determine all the 



coefficients of eq. LI was missing. 

Later it was shown by Britto Cachazo and Feng || that the dj coefficients can be 
easily separated from the rest and computed by introducing quadruple cuts in which the 
loop integration momentum is completely frozen by four on-shell conditions. It was then 
possible to perform a full reconstruction of the amplitude for theories with only boxes, 
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such as N = 4 super- Yang-Mills. However, a systematic procedure to get all the other 
coefficients was still missing. 

Recently the problem of determining in a systematic way the coefficients di, Ci, bi and 
Oj was completely solved by the OPP method of refs. || and [10|. Within this method, 



eq. [LT| is substituted by its unintegrated counterpart, at the price of introducing the so 
called spurious terms, defined by the property of vanishing upon integration over the loop 
momentum q. In practice, since the functional form in q of the spurious terms is universal, 
one has to find, besides di, Ci, bi, and aj, an additional set of coefficients. The OPP method 
allows to find all those coefficients by computing the unintegrated amplitude at different 
values of q for which 4, 3, 2 and 1 propagators vanish. At each stage, the coefficient that 
have been already computed are numerically subtracted from the original amplitude, so 
that, by using such an OPP subtraction, it is possible to disentangle all the coefficients in 
a systematic way. The OPP approach was inspired by the unitarity method and the tensor 
reduction at the integrand level [pd|]. 

More recently the OPP subtraction method was used by the authors of [O] together 
with the Unitarity approach, giving rise to the so called generalized Unitarity techniques, 
that, nowadays, include both semi- analytical |l3|] and fully numerical versions |14|- p5| . 
Nevertheless, in practice, only the so called cut-constructible part of the amplitude, namely 
that one proportional to the one-loop scalar functions, can be easily obtained. The remain- 



ing Rational Terms R |lq1 -|17l require some additional work. For instance, in [[[J] the Ra- 
tional part is obtained by explicitly computing the amplitude at different integer values of 
the space-time dimensions (n). Other possibilities are to get them through n-dimensional 
cuts [15] or with the help of recursion relations fi~9[| . 

On the other hand, in the OPP method, two classes of terms contributing to R nat- 
urally arise [20]. The first class, called Ri, can be derived straightforwardly within the 



same framework used to determine all other coefficients, while the second class, called R2, 
is coming from the (n — 4)-dimensional part of the amplitude and can be obtained by 
computing, once for all, tree-level like Feynman Rules for the theory under study. More- 
over, it is worthwhile to mention that only the full R = Ri + R2 constitutes a physical 
gauge-invariant quantity in dimensional regularization. On the other hand, Ri can be di- 
rectly read out from the analytic expressions of the cut-constructible part of the amplitude, 
irrespectively of the method used to derive it. 

In this paper, we explicitly compute the entire set of Feynman Rules producing R2 
needed in any (massive or massless) QCD 1-loop calculation. We perform our calculation 
in the £ = 1 't Hooft-Feynman gauge. As a consequence, also the 5Z couterterms [Q] needed 
to build renormalized scattering amplitudes should be computed in the £ = 1 gauge. 

In the next section we briefly recall the origin of R2 and give a detailed computational 
example. In section ||, we list our results and present numerical comparisons with known 
amplitudes. In section || we draw our conclusions and, in three appendices, we collect 
diagrams and formulae used for the calculation as well as our results expressed in the color 
connection language. 



- 2 - 



2. The origin of R 2 



Before carrying out our program, we spend a few more words on the origin of R2, that is 
also necessary for setting up the framework of our calculation. 

Our starting point is the general expression for the integrand of a generic m-point 
one-loop (sub-)amplitude 

A(q) = N{q) , A = (? + ft) 2 -"i?, (2-1) 

L>oL>i • • • Mn-1 

where q is the integration momentum. In the previous equation, dimensional regulariza- 
tion is assumed, so that we use a bar to denote objects living in n = 4 + e dimensions 
and a tilde to represent e-dimensional quantities. Notice that, when a n-dimensional in- 
dex is contracted with a 4-dimensional (observable) vector u„, the 4-dimensional part is 
automatically selected. For example 

q • v = (q + q) • v = q ■ v and f = 7^ = f . (2-2) 

An important consequence is 

t = q 2 + f • (2.3) 

The numerator function N(q) can be further split into a 4-dimensional plus an e-dimensional 
part 

N(q)=N(q)+N(q 2 ,q,e). (2.4) 

N(q) lives in 4-dimensions while N(q 2 ,q,e), once integrated, gives rise to the RTs of kind 
R2, defined as 

R2= (^y d q D D 1 ...D m _ 1 - (2 - 5) 

To investigate the explicit form of N(q 2 ,q,e) it is important to understand better the 
separation in eq. |2.4| . From a given integrand A(q) this is obtained by splitting, in the 
numerator function, the n-dimensional integration momentum q, the n-dimensional gamma 
matrices 7^ and the n-dimensional metric tensor g^ y into a 4-dimensional component plus 
remaining pieces: 

Q = Q + Q, 

5 AP = <r+£r. (2.6) 

A practical way to determine R2 is then computing, once for all and with the help of 
eq. ^1], tree-level like Feynman Rules by calculating the R2 part coming from one-particle 
irreducible amplitudes up to four external legs. The fact that four external legs are enough 



is guaranteed by the ultraviolet nature of the RTs, proven in [16]. Through eq. 2.5 a set 
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Figure 1: Diagrams contributing to the gluon self-energy. 

of basic integrals with up to 4 denominators is generated, containing powers of q and e in 
the numerator. A list that exhausts all possibilities in the £ = 1 't Hooft-Feynman gauge 
is presented in appendix [A|. Notice that, according to the chosen regularization scheme, 
results may differ. In eq. |2.5| we use the 't Hooft-Veltman (HV) scheme, while in the Four 
Dimensional Helicity scheme (FDH), any explicit e dependence in the numerator function 
is discarded before integration. Therefore 



R 3 
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FDH (2ttY 



N(q 



r.2 



0) 



(2.7) 



m— 1 



As an explicit and simple example of the described procedure, we detail the calculation 
of R2 coming from the gluon self-energy. The contributing diagrams 1 are drawn in fig. [I]. 
As for the ghost loop with 2 external gluons, we can write the numerator as 



N(q) 



9' 



f aibc f a2ch {p + q)^q^ . 



2.8) 



(2vr) 4 

Since [i\ and 112 are external Lorentz indices, that are eventually contracted with 4 



dimensional external currents, their e-dimensional component is killed due to eq. 2.2 
Therefore, R2 = for this diagram, being N(q 2 ,q,e) = 0. With this same reasoning, 
one easily shows that ghost loops never contribute to R2, even with 3 or 4 external gluons. 

The contribution due to Nf quark loops is given by the second diagram of fig. |l|, whose 
numerator reads 



N(q) 



9 



(M 



-N 



(2.9) 



where the external indices \i\ and \i2 have been directly taken in 4 dimensions. By anti- 
commuting 7^ 2 and $ and using the fact that, due to Lorentz invariance, odd powers of q 
do no contribute, one immediately arrives at the result 



N(q 2 ) = ^4 N f 5 aia2 g^ 1U2 q 2 



(2.10) 



Eq. 2.10, integrated with the help of the first one of eqs. A.l, gives the term proportional 



to Nf in the 2-point effective vertex of fig. |2|. 



1 0ur conventions and notations are listed in Appendix [b| 
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Figure 2: Effective vertices contributing to R2 in pure QCD. X)p(234) stands for a summation over 
the six permutations of the indices 2, 3 and 4, and {t ai t a i} = t ai t a i + t a H ai . Xhv = I in the HV 
scheme and \hv — in the FDH scheme. N co i is the number of colors and Nf is the number of 
fermions running in the quark loop. 
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Figure 3: Effective vertices contributing to R2 mixed QCD. Xhv — 1 in the HV scheme and 
Xhv = in the FDH scheme. In the case of neutral external vectors or scalars, the formulae should 
be read as the contribution given by one quark loop. In the case of charged external particles, they 
refer instead to the contribution of one quark family. 
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Finally, the numerator function of the diagram with a gluonic loop reads 

N{q) = -^^f aibC f a2C X^(P-^-P^)V^(-p,-q,q + P ), (2.11) 
with V given in eq. [B.l| . The contraction of the two V tensors gives terms containing 

q 2 = q 2 +f, and g &p g^ = A + e, (2.12) 

from which one obtains 

Q 2 r 

N(q 2 ,q,e) = -^"f'&f 1 ** [ 2^ 1W § 2 + 4c Wm2 + 2eg w p w + 2cp w ? w + ep w p w 

(2-13) 

Performing the integration, gives the expression written in the first line of fig. 0, where the 



contribution proportional to Xhv is generated by the e dependence of eq. |2. 13 . 

The complete set of effective vertices obtained with the described technique is presented 
in the next section. 

3. Results 

As already explained, 1-loop irreducible Feynman diagrams up to 4 external legs are suffi- 
cient to compute R2 for any amplitude with any number of external legs. Each contributing 
diagram has been calculated analytically by using the Feynman rules listed in appendix [B], 
which also contains the list of the relevant graphs. Different contributions have then been 
summed and reorganized to identify the effective R2 vertices listed in figs. ^ and |3], which 
represent the main result of this work and that allow to determine R2 needed in the com- 
putation of the NLO QCD corrections to any process in the Standard Model. 

In fig. |2] we collect the "pure" QCD effective vertices, namely all vertices generated 
by QCD corrections to processes with external QCD particles. The complete set of con- 
tributing diagrams is given in fig. [|. In fig. ^ we list, instead, the "mixed" QCD vertices 
generated by QCD corrections to processes containing at least one external electroweak 
particle. In this paper, this last class is parametrized by introducing generic couplings of 
a (pseudo)-vector and of a (pseudo)-scalar with a quark line, as in the last two vertices of 
fig. [|. The non vanishing diagrams contributing to R2 are listed in fig. ||. 

In all figures, N co i is the number of colours, Nf is the number of fermions running 
in the quark loop and Xhv is a parameter allowing to read our formulae in two different 



regularization schemes: Xhv = 1 corresponds to the HV scheme of eq. 2.5 while Xhv = 
in the FDH scheme defined in eq. |2.7| . 

Notice that the whole structure of the three-gluon effective vertex is always propor- 
tional to the tree-level, while the four-gluon effective vertex is more complicated. 

Notice also that, when a completely antisymmetric e tensor occur in the formulae of 
the mixed QCD vertices, it always multiplies the axial coupling a. Therefore, summing 
over all quark loops gives zero in the Standard Model, due to the anomaly cancellation. 
Such terms can then be taken to be zero from the very beginning. 
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In figs. |2] and || the effective vertices are given in terms of traces of color matrices 
and structure constants. The same result in terms of color connections is presented in 
appendix ^. 

Just as a showcase of our ability to reproduce the rational terms R2 correctly for 
higher multiplicity of external legs, we have computed the 1-loop six gluon amplitude using 
an extension |21] of HELAC-PHEGAS [22|, HELAC-lloop, that includes virtual corrections 



through an interface with CutTools lCj. 



The comparison against the results of ref. [23] is given in table [T], that refers to one 



color configuration only and to the following phase space point: 

pi = (-3.000000000000000, 1.837117307087384, -2.121320343559642, 1.060660171779821) 

p 2 = (-3.000000000000000, -1.837117307087384, 2.121320343559642, -1.060660171779821) 

p 3 = (2.000000000000000, 0.000000000000000, -2.000000000000000, 0.000000000000000) 

Pa = (0.857142857142857, 0.000000000000000, 0.315789473684211, 0.796850604480708) 

p 5 = (1.000000000000000, 0.866025403784439, 0.184210526315789, 0.464829519280413) 

p 6 = (2.142857142857143, -0.866025403784439, 1.500000000000000, -1.261680123761121) (3.1) 

We find an excellent agreement among the results when including the R2 contribution. 
Finally, we mention that, always with the help of HELAC-lloop, we successfully compared 
our predictions for six-quark 1-loop amplitudes (with three different flavours) with the 
results produced by the GOLEM group [fl. In addition, for all sub-processes included in 



the 2007 Les Houches wish list |24], we explicitly checked that the validity of the Ward 



Identity for a single on-shell external gluon (when present) is preserved by the sum Ri + R2. 



4. Conclusions 

We have derived the tree-level Feynman rules needed to compute the Rational Terms R2 
in QCD, both using explicit color configurations and in the color connection language. We 
listed all effective vertices generated by QCD corrections to processes with external QCD 
particles and all possible mixed QCD effective vertices generated by QCD corrections to 
processes with at least one external EW particle. The inclusion of the derived vertices in 
an actual calculation gives numerical agreement with known expressions for processes up 
to 6 external legs. So we have explicitly checked that 1,2,3 and 4-point vertices are enough 
to solve the problem for an arbitrary number of external legs. In addition, all relevant 
integrals needed to compute R2 in the £ = 1 't Hooft-Feynman gauge have been explicitly 
listed. The next obvious step is the determination of the Feynman Rules needed in the 
complete Standard Model. We leave this for a future publication. 
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Table 1: Results for the finite part of the 1-loop virtual amplitudes for some helicity configurations 
for the case of six external gluons for the phase space point given in the text. The first row for each 
helicity configuration is the tree-order result. The second (unitarity) and the third (semi-numerical) 
rows are the results for \A$\ taken from [^3|. The fourth is our result for the cut-constructible, CC 
(with renormalization scale /i = v^); ^i an d R2 terms in the HV scheme as well as for the \A§\ 
in FDH scheme, to facilitate comparisons. The relation of HV scheme result to the FDH scheme 
result is given in [^3| . 
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Calculations for Physics Discoveries at Colliders). The research of R.P., M.V.G. and P.D. 
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Appendices 

A. The needed integrals 

In this appendix, we collect the integrals needed to perform our calculation of R2. A non 
vanishing contribution is generated only by integrals with zero or higher superficial degree 
of divergence. They fall in two classes, namely integrals involving even powers of q (odd 
powers do no contribute due to Lorentz invariance) and Pole Parts (P.P.) of ultraviolet 
divergent integrals. This second class is relevant when using regularization schemes (such 
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as the HV one) where the e dependence in the numerator is kept. In the following, we 
further classify the integrals according to the number of denominators Di = (q+Pi 
The results for the Pole Parts have been checked against ref. 0]. 

2-point integrals: 
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3-point integrals: 
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with p ijk =pi+ Pj +p k . 
4-point integrals: 
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B. QCD Feyman Rules and diagrams 

In this appendix, we present the Feynman rules and the diagrams used in the calculation. 
In fig. El we list QCD propagators and vertices as well as our parametrization of the Vqq 
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and Sqq couplings. Ghosts are drawn with dashed arrows, vectors with wavy lines and 
scalars with dotted lines; Greek letters denote Lorentz indices; k, I = 1,2,3 are the three 
colors of the quarks while all remaining color indices range from 1 to 8; f abc is the QCD 
SU(N co i) structure constant and t a (a = 1, • • • , 8) are the color matrices in the fundamental 
representation; m q is the quark mass and V Ulfl2fl3 (pi,p2,P3) is given by 

V I mi 1 *2 I m 3 (P1,P2,P3) ^/n^te-PlU +6>2/«(P3 -P2) W + Wife -P3)/« • (B- 1 ) 

Finally, in fig. [5] and |(| we draw the pure QCD graphs and the mixed QCD diagrams which 
give a non vanishing contribution to R2. As explained in section 0, diagrams involving QCD 
FP ghosts do not to contribute to R2 and are not included. 
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Figure 4: Feynman rules used for the computation. The last two vertices parametrize a generic 
coupling of a (pseudo) -vector V and of a (pseudo)-scalar S with a quark line, respectively. 
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Figure 5: Pure QCD diagrams used for the calculation. 



C. Effective vertices in the color connection language 



In this appendix, we write down the Feynman rules for the QCD effective vertices in the 
color connection language. Such rules are obtained by contracting any gluon index dj, 
appearing in the vertices of figs. |2| and [3| , by a color matrix t^.. Any gluonic color index 
a,i is therefore projected out in terms of two quark like color and anti-color indices fcj and 
li. By then summing over gluon indices with the rule 



j.a j.a 
z kl l ij 



$kj5il 



1 



Ncol 



(C.l) 



the color part of the effective vertices can be entirely written down in terms of <5's, which 
correspond to color connections. Graphically, a color connection can be represented with 
a solid line, in such a way that two solid lines stand for a gluon, while one single solid 
line symbolize a quark. Finally, different color lines can be connected by the exchange of 
a scalar colorless gluon, represented by a dashed line. In such a language, the pure QCD 
effective vertices of fig. |2|can be written as in figs. [?] and ||. Analogously, the last five mixed 
QCD vertices of fig. || give the results reported in figs. ^-11. 
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Figure 6: Mixed QCD diagrams contributing to R2. 
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